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Abstract. I summarize four of the most important areas of uncertainty in the study of the chemistry and cooling of gas with 
zero or very low metallicity. These are: i) the importance and effects of HD cooling in primordial gas; ii) the importance of 
metal-line and dust cooling in low metallicity gas; iii) the impact of the large uncertainties that exist in the rate coefficients of 
several key reactions involved in the formation of H2; and iv) the effectiveness of grain surface chemistry at high redshifts. 
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INTRODUCTION 

Astrochemistry - chemistry in an astrophysical context 
- is a huge subject, owing to the extremely large num- 
ber of reactions that are possible, and the complexity of 
the resultant reaction networks. For instance, the reaction 
sets used to study the chemistry of local molecular clouds 
typically include more than 300 reactants and more than 
4000 reactions simply to represent the gas-phase chem- 
istry. The inclusion of grain surface chemistry makes 
things even more complex. However, when studying the 
formation of the earliest generations of stars, our primary 
concern is to properly model the thermal behaviour of the 
gas, as this strongly influences its hydrodynamical evo- 
lution. We can therefore focus our attention solely on the 
chemistry of the major coolants, allowing us to sidestep 
much of the chemical complexity, and to include only 
those reactants and reactions in our model that are neces- 
sary for understanding the thermal evolution of the gas. 

An influential study of the chemistry of metal-free gas 
performed along these lines by [ 1 ] demonstrated that out 
of the many hundreds of possible reactions, only 28 were 
needed to accurately model the thermal evolution of the 
gas over a wide range of parameter space. Subsequently, 
other authors have added a few additional reactions to 
this simplified model that are very important in some 
circumstances (e.g. HD formation and destruction, J2l; 
three-body H2 formation, Q), but not in others. 

A similar approach has recently been applied to metal- 
poor gas by [4] and [5], under the assumption that the 
dominant coolants in metal-poor gas are the same as 
those in the local atomic and molecular ISM (i.e. carbon, 
oxygen and silicon fine-structure lines, CO and H2O 
rotational and vibrational emission, and dust). In this 
case, the required number of reactions is an order of 
magnitude larger than in the metal-free case, but is still 
small enough to be computationally tractable. 

Nevertheless, despite our success at identifying the 



most important reactions and hence simplifying the treat- 
ment of the chemistry to the point where it can be mod- 
elled self-consistently within large, three-dimensional 
hydrodynamical simulations, a number of issues remain 
unresolved. In this contribution, I briefly discuss a few of 
what I consider to be the most important areas of ongo- 
ing uncertainty or confusion in the study of the chemistry 
and cooling of metal-free and metal-poor gas. 



PROBLEM 1: THE ROLE OF HD 
COOLING 

At temperatures 200 < T < 10000 K in a metal-free 
gas (or in low-density metal-poor gas; |01), molecular 
hydrogen is the most abundant molecular species and is 
also the dominant coolant. However, at low temperatures, 
H2 cooling becomes increasingly ineffective, owing to 
the large energy difference that exists between its / = 
and 7 = 2 rotational levels, and the fact that the lower 
energy / = 1 — > transition is strongly forbidden. At 
these low temperatures, HD cooling is far more effective, 
in the sense that the cooling rate per molecule of HD is 
much larger than the cooling rate per molecule of H2. 
However, the small size of the cosmic D/H ratio means 
that in general HD is far less abundant than H2. Because 
of this, for a long time it was unclear whether or not HD 
was ever an important coolant in practice. In the past few 
years, however, considerable effort has been devoted to 
exploring the consequences of HD cooling in the early 
Universe, and its role has now become much clearer. 

In H, the impact of HD cooling in the first star- 
forming protogalaxies was examined, and it was shown 
that given a D/H ratio consistent with observationally 
measured values, HD cooling was unimportant. More 
recently, the role of HD cooling within a wider range of 
protogalaxies was examined ([8]), with the conclusion 
that in very low mass protogalaxies (M < 3 x 10 5 M Q ), 



HD cooling may be marginally important. 

These calculations typically adopted an initial tem- 
perature and chemical composition for the gas consis- 
tent with the values applicable to the undisturbed IGM. 
However, other studies that start with initially ionized gas 
(e.g. dHljli) find that in this case, HD cooling is impor- 
tant, and that it can often cool the gas all the way down 
to the temperature floor set by the CMB. 

Why do these two sets of calculations differ so much 
on the importance of HD? The answer lies in the low 
temperature chemistry of the HD molecule. In most cir- 
cumstances, its formation and destruction are dominated 
by the ion-neutral exchange reactions: 



H 2 +D H 
HD + H H 



HD + H H 
H 2 + D+ 



(1) 

(2) 



Reaction [TJ is exothermic, but reaction [2] is endothermic 
by 462 K (lHHl). As a result, at low gas temperatures, 
considerable chemical fractionation occurs: the equilib- 
rium value of the HD/H2 ratio is enhanced over the cos- 
mic D/H ratio by a factor ~ e 462 / r . The effect of this on 
the cooling of the gas is illustrated in Figure [TJ where 
I plot the cooling rate of the gas as a function of tem- 
perature for gas with n = 100 cm~ 3 and jc(H 2 ) — 10~ 3 , 
and where x(HD)/x(H 2 ) is calculated assuming that re- 
actions [TJ and 12 have reached equilibrium. In practice, 
fractionation may not proceed fully to equilibrium ow- 
ing to the depletion of the necessary D + ions from the 
gas by recombination, so this plot should be regarded as 
giving an upper limit on the effectiveness of HD cooling 
at this density. Figure QJ demonstrates that at T ~ 150 K, 
the increase in the HD abundance due to fractionation 
more than offsets the decrease in the H2 and HD cooling 
rates per molecule. Consequently, gas that can cool to 
this critical temperature will begin to cool even more ef- 
fectively, and thus can easily reach temperatures as low 
as Tcmb- The critical temperature that must be reached 
for HD cooling to take over in this fashion varies with n, 
but is never larger than 250 K, and exceeds 200 K only 
for gas with n > 10 4 cm~ 3 . 

It is now straightforward to explain why HD cool- 
ing becomes important in gas cooling from an initially 
highly ionized state, but not for gas cooling from a low- 
ionization state. In the former case, it has long been 
recognised that the out-of-equilibrium nature of hydro- 
gen recombination in the cooling gas leads to enhanced 
production of H 2 . As a result, gas cooling from a high- 
ionization state can reach a lower temperature than gas 
cooling from a low-ionization state. Although the tem- 
perature difference is not large, it turns out that in the for- 
mer case, the gas can usually cool to T alX , at which point 
HD takes over, while in the latter case it falls somewhat 
short, typically never cooling below T ~ 200 K. 
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FIGURE 1. Cooling rates of H 2 (solid line) and HD (dashed 
line) in a gas with atomic hydrogen number density «h = 
100 cm~ 3 . The fractional H 2 abundance is fixed at x(H 2 ) = 
10~ 3 , while the HD/H 2 ratio, x(HD)/x(H 2 ), is calculated as- 
suming that reactions [TJ and [2] have reached equilibrium. The 
maximum abundance of HD is constrained to be no larger than 
the total deuterium abundance, here taken to be x(D,tot) = 
2.6 x 10~ 5 . This constraint is responsible for the sudden down- 
turn in the HD cooling rate at low T . 



This explanation accounts for the differences in be- 
haviour observed in previous numerical simulations with 
different initial conditions, but also prompts a couple 
of further questions. First, what effect does the pres- 
ence of an ultraviolet background in the Lyman- Werner 
bands of H 2 have on the ability of the gas to cool be- 
low r cr it? Most previous simulations that include the 
effects of HD cooling have assumed that the ultravio- 
let background is zero, or at least is too small to be 
significant. Recently, however, the impact of a Lyman- 
Werner background on the effectiveness of HD cool- 
ing was investigated in 111 211 using a one-zone chemi- 
cal model. This investigation suggests that if the field 



strength is greater than3 x 10 _zz ergs _1 cm _2 Hz _1 sr _1 , 
then the gas never cools below 200 K, and so HD cooling 
never becomes important. However, this calculationsuf- 
fers from the problem - common to all one-zone mod- 
els - that the dynamical evolution is not solved for self- 
consistently with the thermal evolution. Further investi- 
gation of this issue using a proper hydrodynamical ap- 
proach would be extremely useful. 

Second, does the extra cooling at low temperatures 
provided by the HD actually promote fragmentation? 
Judging by recent work (e.g. [10]), the answer is no: 
the additional cooling reduces the fragment mass scale, 
but despite this only a single fragment forms in each 
halo. However, the sharp rise in the cooling rate at T < 



?crit apparent in Figure Q] suggests that HD cooling may 
render the gas thermally unstable in some circumstances, 
which would tend to promote fragmentation. 

PROBLEM 2: THE IMPORTANCE OF 
METALS 

One of the most important unsolved problems in the 
study of star formation in very low metallicity gas is the 
role that cooling from gas-phase metals and dust grains 
plays in determining the stellar IMF. A popular hypothe- 
sis (introduced in IU3I0 holds that it is the enrichment of 
the star-forming gas with metals to a level greater than 
a critical metallicity Z cr j t that brings about a transforma- 
tion in the IMF from one dominated by massive stars to 
one more closely resembling the standard Salpeter IMF. 
The causal link between metal enrichment and the for- 
mation of lower mass stars is the greater effectiveness of 
metals as coolants (particularly at low temperature) com- 
pared to H2, which is supposed to enable the gas to grav- 
itationally fragment down to much smaller mass scales. 

Two different forms of metal-related cooling have 
attracted considerable attention: fine-structure cooling 
from C and O (and in some cases also Si and Fe; [14]), 
which is effective primarily at densities n < 10 6 cm~ 3 
dill EH, and dust cooling, which is effective at 
much higher densities lfl7l ITsl [l9[|20ll . 

The effectiveness of fine-structure cooling at promot- 
ing fragmentation remains highly uncertain. The original 
simulations of IU3I1 found fragmentation at Z = 1O~ 3 Z 
but not at Z = 1CT 4 Z , implying that 1CT 4 < Z crit < 
1O~ 3 Z . However, these simulations did not include the 
effects of H2 cooling. Very similar simulations have re- 
cently been performed that do include H2 cooling ( J2lll ). 
In this case, fragmentation occurs even in runs with Z 
0. suggesting that the fragmentation seen in this type of 
simulation owes more to the choice of initial conditions 
than to metal-enrichment. Furthermore, simulations per- 
formed using alternative initial conditions in which the 
gas is hot and ionized ( 112 ill ) find no fragmentation at 
Z < 1O~ 3 Z , whereas simulations by [16], performed 
using yet another set of initial conditions, do find frag- 
mentation at Z = 10~ 3 Z . The logical conclusion to be 
drawn from these contradictory results is that whether 
or not metal-enriched gas fragments depends to a large 
extent on its initial state, and that without a better under- 
standing of this, one cannot be certain about the impor- 
tance of fine-structure cooling in driving fragmentation. 

In the case of dust cooling, recent hydrodynamical 
simulations (described in lfl9fl . l2(ill and elsewhere in 
these proceedings) find that the gas readily fragments 
following the onset of dust cooling, provided that the 
available dust surface area per unit gas mass exceeds one 



part in 10 5 of the value in solar metallicity gas. Although 
the simulations of lfl9ll . which consider the special case 
of gas with zero angular momentum, find only limited 
fragmentation, the simulations of [20], which consider 
the more general case of non-zero angular momentum, 
find a high level of fragmentation, with ^100 distinct 
fragments formed within only a few hundred years of the 
formation of the first object. 

Nevertheless, despite these recent results, a number of 
serious questions remain. To begin with, both sets of sim- 
ulations were performed using a tabulated equation of 
state, based on the thermal evolution of the gas found 
in one-zone models of low-metallicity protostellar col- 
lapse (e.g. [18]). However, one-zone models of this type 
do not self-consistently couple the thermal and dynam- 
ical evolution of the gas, and also cannot properly ac- 
count for effects such as rotation. Conclusions based on 
one-zone models of collapse in metal-free gas have in the 
past proved to be highly misleading (see e.g. the discus- 
sion in [22]). Moreover, by using a tabulated equation of 
state, we implicitly assume that the gas temperature ad- 
justs instantaneously to any changes in density; in prac- 
tice, such changes are not instantaneous, but instead take 
place over a thermal timescale. It is possible that by mak- 
ing this implicit assumption, we are overestimating the 
propensity of the gas to fragment. Work aimed at clarify- 
ing these issues by solving the full thermal energy equa- 
tion self-consistently within the hydrodynamical simula- 
tions is currently ongoing. 

The other major concern about the current results is 
that they do not take into account feedback in the form 
of radiation from the accretion shocks surrounding the 
newly formed protostars. In local star-forming regions, 
this form of feedback may play a crucial role in stabiliz- 
ing massive cores against sub-fragmentation ( 12311 ). and 
so its effects should clearly also be included in the high 
redshift simulations. However, there are reasons to be- 
lieve that its effects will be far less significant in a high- 
redshift, low-metallicity context. For one thing, the gas 
in these low-metallicity cores that cools rapidly and frag- 
ments is typically not optically thick, and so will absorb 
only a fraction of the radiation emitted by the accret- 
ing protostars. For another, the temperature of the dust 
in these objects is very much higher than in local star- 
forming regions (Td ~ 200 K at the relevant densities, 
compared to T& ~ 10 K locally), and so a much higher 
energy input is required to significantly alter the temper- 
ature of the gas. For these reasons, it seems unlikely that 
the inclusion of feedback will significantly alter the out- 
come of the simulations, but this expectation must nev- 
ertheless still be properly tested. 



PROBLEM 3: RATE COEFFICIENT 
UNCERTAINTIES 

Our ability to construct accurate models of the chemical 
evolution of metal-free or metal-poor gas is constrained 
by the level of accuracy to which the rate coefficients of 
the key chemical reactions have been determined. Most 
of the important reactions involved in the chemistry of 
H9 and HD have well-determined rate coefficients (JTt], 
12411 V However, there are a few important exceptions. 

Two of these exceptions are the associative detach- 
ment of H with H 

H +H^H 2 + e , (3) 

and the mutual neutralization of H with H + 

H +H+^H + H. (4) 

The available experimental and theoretical data on the 
low temperature (T < 10 4 K) behaviour of these two 
reactions is summarized in B25I1 . In both cases, the rate 
coefficients are uncertain by at least a factor of a few. 
In gas with a low fractional ionization, this uncertainty 
is unimportant, as reaction [3] occurs at a much faster 
rate than reaction regardless of which of the possible 
values are selected for the two rate coefficients. In gas 
recombining from an initially ionized state, on the other 
hand, the uncertainties in reactions [3] and [4] combine to 
render the H2 abundance uncertain by as much as an 
order of magnitude at some points during the chemical 
evolution of the gas. Some possible consequences of this 
were also explored in 1I25I1 . 

Another example is the large uncertainty that exists 
in the rate coefficient of the three-body H2 formation 
reaction 

H + H + H^H 2 + H, (5) 

the dominant source of H 2 in gas at densities n > 
10 8 cm~ 3 . In Figure |2j I show a selection of rate co- 
efficients quoted in the literature for this reaction. The 
corresponding numerical values and references are listed 
in Table [1] This compilation includes one new rate co- 
efficient (F), presented here for the first time. This was 
computed using the same detailed balance approach as in 
l26ll . but illustrates the effect of using an H 2 collisional 
dissociation rate taken from l27ll in place of the rate from 
JH used in JH]. 

It is clear from Figure [2] that although there is reason- 
able agreement between many of the rates at 5000 K, 
there is a substantial uncertainty in the rate in the temper- 
ature range 200 < T < 2000 K relevant for population III 
star formation. Moreover, there is no sign that the rates 
are converging: two of the most recent determinations (A 
and D) are the two most widely discrepant rates. 



CM 

I ______ 

o I 




2 100 1000 10 4 

Temperature (K) 

FIGURE 2. Three-body H 2 formation rate coefficients, plot- 
ted as a function of temperature. Key: A - lower solid line; B 
- dashed line; C - dot-dashed line; D - dotted line; E - star 
symbols; F - upper solid line. 



TABLE 1. Rate coefficients for three-body H 2 formation 



ID 


Rate coefficient (cm 6 s 1 ) 


Reference 


A 


1.14 x io- 31 r-°- 38 


T < 300K 


[_ 




3.9 x io- 30 r- 10 


T > 300K 




B 


5.5 x 




[29J 


C 


8.8 x 1(T 33 




[30] 


D 


1.44 x icr 26 r- L54 




[26] 


E 


1.4 x 1(T 31 


T = 3000 K 


[31J 




8.2 x 1(T 32 


r = 5000K 


F 


7.7 x i(T 31 r- a464 




This work 



The 7i 3 density dependence of the three-body H 2 for- 
mation process ensures that in collapsing primordial pro- 
tostellar cores, the gas will rapidly become fully molec- 
ular at high densities, regardless of the rate coefficient 
adopted. However, the choice of rate coefficient will 
strongly influence the three-body H 2 formation heating 
rate, which is a major energy input into the gas at these 
densities, and may therefore also affect the further dy- 
namical evolution of the gas. Simulations aimed at ex- 
ploring the consequences of this uncertainty are currently 
in progress. 

PROBLEM 4: GRAIN SURFACE 
CHEMISTRY 

The key role that dust grains play in the chemistry of 
molecular gas in the Milky Way and in other local galax- 
ies has been understood for more than 40 years. Very few 



of the H2 molecules in these systems were formed in the 
gas phase via the H or pathways that are so im- 
portant in primordial gas. Instead, most were formed by 
reactions between hydrogen atoms adsorbed on the sur- 
face of interstellar dust grains. Moreover, H2 is far from 
the only molecule that can form in this manner. 

However, the importance of grain surface chemistry 
at high redshift remains an open question. There are 
two key issues. The first is that the amount of dust 
present in high redshift, low metallicity systems is not 
well constrained. A number of calculations have been 
performed that explore the formation of dust grains in 
high-redshift supernova remnants produced by metal- 
free progenitors ((H SHil Q), but the results 
of these calculations do not agree on the details of the 
resulting grain size distribution, and so the available 
grain surface area per unit gas mass remains uncertain. 

The second key issue is the efficiency of molecule for- 
mation on high redshift dust grains. On local dust grains, 
most chemistry occurs through the interaction of van der 
Waals bonded (or 'physisorbed') atoms. The evaporation 
timescale for these atoms is typically orders of magni- 
tude longer than the accretion timescale, and so the atoms 
have plenty of time in which to react. At high z, how- 
ever, the dust grains are warmer, owing to the effects of 
heating by the CMB, and so the evaporation timescales, 
which depend exponentially on Td ust , are vastly smaller. 
For important atoms such as H, C or O, this means that 
tevap <S face, and so the likelihood of there being more 
than a single atom physisorbed on a given grain at any 
given time is extremely small. Taken at face value, this 
implies that the efficiency of molecule formation on high 
redshift dust grains should be extremely small. 

If one applies a similar argument to warm dust in lo- 
cal photodissociation regions, then one predicts a sim- 
ilar outcome: molecule formation, and in particular the 
formation of H2, should be extremely inefficient there. 
However, this conclusion conflicts with observationally- 
derived constraints on the H2 formation rate in these re- 
gions, suggesting that some or all of the H2 formed in 
these regions is formed from H atoms that are chem- 
ically bonded (or 'chemisorbed') to the surface of the 
grains ([37]). A detailed model for the formation of H2 
on grains that have both physisorption and chemisorp- 
tion binding sites has recently been constructed by IB8I1 
and applied to the question of high-redshift H2 forma- 
tion by [39]. The results of this model suggest that once 
chemisorption is taken into account, H2 formation re- 
mains efficient even at very high redshifts. However, this 
model is not without its critics (e.g. |40|]). and as yet no 
similar modelling effort has been completed for the for- 
mation of any species other than H2, although work along 
these lines is currently under way. 
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